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Tripod-scheme cold atoms interacting with laser beams have attracted considerable interest for 
' their role in synthesizing effective non-Abelian vector potentials. Such effective vector potentials can 

be exploited to realize an all-optical imprinting of geometric phases onto matter waves. By working 
5— ( , on carefully designed extensions of our previous work, we show that coherent lattice structure of 

cold-atom sub-wavepackets can be formed and that the non-Abelian Aharonov-Bohm effect can be 
' easily manifested via the translational motion of cold atoms. We also show that by changing the 

frame of reference, effects due to a non-Abelian vector potential may be connected with a simple 
dynamical phase effect, and that under certain conditions it can be understood as an Abelian 
geometric phase in a different frame of reference. Results should help design better schemes for the 
control of cold-atom matter waves. 
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c3 ■ I- INTRODUCTION 

^3 . 

^ / As one of many significant developments in using cold-atoms to achieve quantum simulations, tripod-scheme cold 
atoms interacting with three laser beams have attracted considerable interest for their role in synthesizing effective 
CsJ non-Abelian vector potentials [ll-@|- Specifically, for cold-atoms with a tripod internal level structure and in the 
^ presence of plane-wave laser fields, a simple space-independent (as long as the atoms are well inside all the laser 
. beams) non-Abelian vector potential can be generated. This result has led to a number of interesting predictions, 
' such as negative anomalous reflection Q , cold-atom analog of Datta-Das transistor [9] , negative refraction [l^l and 
cold-atom analog of the so-called Zitterbewegung oscillation (Tll - fTsj . 

Recently, we showed that if the laser beams interacting with a tripod-scheme cold-atom are displaced slowly, then it 
is possible to actively imprint a non-Abelian vector potential induced phase [3] onto the matter wave [15| . Because the 
phases thus obtained are not induced by the translational motion of the cold atom itself, but by active manipulation 



of laser-matter interaction, the imprinting of quantum phases onto matter waves results in a new type of coupling 



between internal and translational motions. Indeed, as shown in Ref. 15[, due to the phase imprinting a wavepacket 
. ^ ' may develop many interesting interference patterns and may even be splitted into many equal- weight copies along 
. a straight line. Because by construction the phase imprinting due to the laser-beam displacement is insensitive to 
the details of the laser beams (such as the laser intensity, actual speed of the moving laser beam, etc), the active 
5^ , displacement of laser beams can potentially provide a robust means of controlling matter wave propagation. Further, 
a direct experimental observation of the phases due to a non-Abelian vector potential can now be reduced to the 
observation of the translational motion of cold atoms. 

The purpose of this work is twofold. First, we consider a few more scenarios of laser-beam displacement, thus 
extending and strengthening our early work |15|. In one scenario where the laser-beams are moved along a square, we 
show that it is possible to split a wavepacket into equal-weight copies on a square lattice. We further show that it is 
also possible to form a triangular lattice of cold-atom sub-wavepackets. In another scenario where the laser beams are 
moved along a circle, we show that one may force a cold-atom wavepacket to "dance" along a circle of a different radius 
or redirect the propagation direction of that wavepacket. As seen below, these new scenarios are carefully designed 
for our second and more important purpose, namely, to address two issues of more fundamental interest. The first 
issue is on a simple but intriguing demonstration of non-Abelian Aharonov-Bohm effect via the translational motion 
of cold atoms. The second issue is on how our perspectives of quantum phases may change as we change the frame 
of reference. We shall explicate that, by changing the frame of reference, i.e., by changing from the laboratory frame 
to a reference frame moving with the laser beams, the effect due to an effective non-Abelian vector potential may be 
connected with a dynamical phase effect, and under certain conditions part of this effect may be understood as an 
Abelian geometric phase obtained from a different frame of reference. These results, based on concrete and explicit 
examples, may enhance our understanding of the dynamics of tripod-scheme cold atoms and motivate experiments 
on the associated control of cold-atom wavepackets. 



2 



This paper is organized as follows. For self-completeness, in Sec. II we provide some necessary details about the 
effective non-Abelian vector potential realized by tripod-scheme atoms. In Sec. Ill we show how moving laser beams 
along a square can lead to formation of a square lattice of cold-atom sub-wavepackets. An extension of this approach 
can lead to the formation of a triangular lattice of cold-atom sub-wavepackets. In analyzing the results we show that 
the physical effects induced by the non-Abelian vector potential can be connected to pure dynamical phases in another 
reference frame. In the same section we also discuss how such laser manipulation can serve as a direct approach to 
the observation of the non-Abelian Aharonov-Bohm effect. In Sec. IV, we provide an interesting relation between an 
evolution matrix as an integral of the non-Abelian vector potential in the laboratory frame and an Abelian geometric 
phase in the frame that moves with the laser beams. In so doing we consider the displacement of laser beams along 
a circle and the resultant motion of cold-atom wavepackets. Section V concludes this work. 



II. NON-ABELIAN VECTOR POTENTIALS REALIZED BY TRIPOD-SCHEME ATOMS 

Tripod-scheme atoms refer to four- level atoms interacting with three laser fields [l6l |. We denote the four internal 
levels as |n), n = — 3. Each of the three transitions |0) -H> |0) O |2), and |0) O |3) is coupled by one laser 
field. This coupling scheme can be realized if states |2), and |3) are degenerate magnetic sub-levels and the three 
coupling fields have different polarizations. For convenience we adopt the same configuration as in Ref. Q, where two 
laser beams are counter-propagating along the a;-axis and the third laser beam is along the z-axis. The associated 
internal Hamiltonian under RWA is given by 



HT^WAA = ^n.n\0){n\+h.C., (1) 

with 



n=l 



= ^h^e-<^^ (2) 
= ^^if^^e^^'-, (3) 

= ^^ocos(Oe"='^ (4) 

where the parameter ^ is set to satisfy cos((^) = V2 — 1, as in Ref. [3], and fc[ is the wavevector of the laser fields. 

The Hamiltonian -ffRWA.4 has two degenerate states with a null eigenvalue. We denote these two degenerate states 
as |£'i(2)) and it is straightforward to find their spatial dependence as follows Q, 

\D^) = - |2))e— '7V2 

\D2) = [cos(0 + |2)) /V2 - sin(0|3)] '^ (5) 

where 

k' = kl[l-cos{0], (6) 

|i) = 1 1) (7) 

|2) = |2)e-*'='-(^-^). (8) 

When the laser field strengths are sufficiently large, can be large compared to any possible two-photon detuning 
induced by fluctuations in laser frequencies and/or Doppler shift. Then if the initial state is in the dark-state subspace 
and if the system parameters are changing slowly, the internal atomic state can evolve within the two-dimensional 
dark-state subspace. As such we focus on the time evolution in the dark-state subspace only. Clearly then, the 
dynamical phase contributed by the internal energy is always zero because i?RWA.4|^i(2)) — 0- Any state therein can 
be spanned as ci\Di) +C2\D2). This expansion henceforth defines a dark state representation. In this representation, 
the mechanical momentum operator becomes 0, Q 



-ihV + h,K{a.j,ex + a^e^), (9) 



where ax.z are Pauli matrices, k — cos(^)fcJ., V represents the gradient in the dark-state representation, and Cx and 
62 are the unit vectors along x and z. As in Ref. Ti], we also introduce an additional constant shift to state |3), which 
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accommodates a detuning of the third laser from resonance by Vs — fi{kl.)^ sin^(^)/2m (m is the mass of atom). Then 
the total effective Hamiltonian becomes 

(■pD\2 

(10) 

This effective Hamiltonian and the explicit form of in Eq. ^ make it clear that a tripod-scheme atom interacting 
with three laser beams effectively synthesize a non-Abelian vector potential, whose a;-component is given by TinUxj 
z-component is given by tiKaz, and y-component is zero. 
The eigenstates of if|f are 

where R is the spatial coordinate (for simplicity the plane- wave normalization factor {27rh)^/^ is not included), and 
\g^) denotes two-component internal-state vectors in the dark subspace. Note that \g^) depends on the direction 
of the wavevector k via its dependence on tpk, the angle between the cc-axis and the atom wavevector k. From Eq. 
([9]), it is straightforward to show that the eigenstates |\E'^'='=) have momenta (k± Kk)h respectively, with their energy 
eigenvalues [{k ± k)K\'^ /2m. Note also that throughout this paper k represents a unit vector along the vector k and 
k represents the modulus of k. 



III. WAVEPACKET CONTROL BY MANIPULATING A NON-ABELIAN VECTOR POTENTIAL 

In Ref. [l^, we proposed to control the matter wave of tripod-scheme atoms by displacing the laser beams along 
a straight line. By considering two different frames of reference, in this section we will first re-analyze our results in 
Ref. We then generalize the scheme in [l^ by displacing the lasers along a square, as illustrated in Fig. 1. As 

another extension we also consider laser displacement along an equilateral triangle. 



A. Displacing laser-beams along a straight line: two perspectives 

1. Perspective from the laboratory frame 

Consider a laser-beam movement along a straight line, e.g., along the path ^4 — i? in Fig. 1. As what we did in Ref. 
[Tsj , it is natural to first analyze this problem in the laboratory frame. Without loss of generality we assume that the 
initial state is given by 

l*!^"^) = |gk(^k = 7r/2))e"'-^ 




(12) 



where we have assumed that the wavevector ko associated with the spatial part of the total wavefunction is zero ■ 
Denoting x = Vdt, where Vd is the speed of the laser beam displacement, then at each spatial point R in the laboratory 
frame the evolution of (ci, C2) is determined by 



k \ ( ci 
n Q )\c2 

^-G.[l), (13) 

where Gx is defined as a 2 x 2 matrix associated with the displacement along the x direction. This expression is 
expected because it refiects the x-component of the non-Abelian vector potential in our dark-state representation 



counter-clockwise 



X 

FIG. 1: One scenario of laser beam displacement along a square for the control of wavepackets of tripod-scheme cold-atoms. 
As explained in the text, a cloud of tripod scheme cold-atoms is simultaneously illuminated by three laser beams. All the laser 
beams are moved either clockwise or counter-clockwise from the starting point A, along the square with the length of each side 
given by a. 



(See Eq. (5)). One tends to call the evolution matrix determined by the above equation as a non-Abelian geometric 
phase. However, because here the involved path in the parameter space is not necessarily closed or cyclic, such an 
evolution matrix is clearly representation-dependent. Nevertheless, if we move the laser beams such that e*'^''-^ = 1, 
then because the parameters for both the internal Hamiltonian and the dark states are subject to cyclic changes, the 
evolution matrix associated with Eq. (13) becomes gauge-independent and therefore can be regarded as a non-Abelian 
geometric phase. In either case, the physical effects discussed below are based on the actual solution to Schrodinger's 
equation and are hence representation-independent. 

In addition to this effect due to the non-Abelian vector potential, a simple dynamical phase also accumulates during 
the laser manipulation. In the laboratory frame, the energy of any state {ci , C2)'^ e^^"^ — {ci\Di) + C2|£'2))e*'"'^ is 
contributed by two terms. The first term is the internal energy of the dark subspace, which is zero at all times. The 
second term is the eigenvalues [(fco ± K)?i]^/2m of the eigenstates given in Eq. ([TT|) . In the case of ko = 0, this term 
is given by 

j^n^^ ^ HMA(2)I4IA(2)))^ ^ (14) 
2m 2rn ' 

which is independent of ci or C2. Hence, so long as the laser displacement is sufhciently slow, the internal state will 
remain in the dark subspace and the kinetic energy will be necessarily given by El^^^ for ko = 0, which is independent 
of ci or 02- 

Based on Eq. (|13p and the dynamical phase determined by i?]^^, one then finds the state evolution during the 
passage A — B, i.e.. 



2 VI/ 2 V-1 



2 V 1 



where 



2 



2m 



is the overall dynamical phase induced by the kinetic energy term E^^^, and 



^2 = VdK 



(17) 
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are two frequencies due to the non-Abelian vector potential. Because, at a time t, the total distance of laser dis- 
placement in the x direction is a; = w^i, the phase factors ±VdKt can also be written as ±kx, depending only on the 
distance of laser displacement. This implies the robustness of this phase acquired by the process, thus manifesting its 
geometric nature. Equation (fT5|) also indicates that the acquired phases zLna at point B are space- independent, and 
hence it will maintain ko = as the wavevector for the spatial part of the total wavefunction. This further justifies 
our treatment here. 

Equation for ko = actually indicates the existence of two group velocities. This can be well understood by 
finding the derivatives of the frequencies a;i(2) with respect to the atomic wavevector fcg in the laboratory frame. To 
that end we consider a nonzero but small Sko in the x direction in the laboratory frame. We then choose another 
frame of reference that moves with atom with a velocity of hSkoix/m. In the new frame of reference the wavevector 
becomes zero again, the laser beams are displaced with a velocity Vd — 'Vd — hSko/m, and Eq. (jlSp can be directly 
applied with the modified moving speed Vd- Since 



dvd _ _ 
dko 

the two group velocities can then be obtained using Eq. ([18 



(18) 



lab d{-Vdn) dvd 

^iG = ^ = — ^ IT = m, 

oko ovd dko 



lab 



duj2 _ d{vdK) dvd 
dko dvd dko 



'fiK/m. (19) 



Because these two group velocities associated with the two components of the evolving state in Eq. (|15|) are different, an 
initial wavepacket of the tripod-scheme cold atoms is expected to split into two, each of the sub-wavepackets possesses 
an internal state "^(Ij 1)'^ or -^(1, —1)^- This result is nicely confirmed by numerical simulation results shown in 

Fig. 2(a). In particular, all our numerical experiments are based on the full Hamiltonian — — h -ffRWA,4 + K , and 
the initial state is chosen as a Gaussian wave packet instead of a plane wave. From Fig. 2(a), it is seen that when 
the point B is reached, two sub-wavepackets located at x = ±!^tB, z = emerge, where tx stands for the time when 
point X is reached. Because the length of square shown in Fig. 1 is given by a, we have ts = a/vd- 

2. Perspective from the frame of reference that moves with the laser beams 

In the frame of reference that moves with the laser beams, called laser frame below, all the three laser beams 
are static by construction, but the atom momentum will be changed. Because the Hamiltonian in Eq. ()10p and 
its eigenstates in Eq. PT|) implicitly assume fixed laser-beams, this laser frame can be treated in a straightforward 
manner. For example, if the laser beams are moving with a speed Vd along the cc-direction, then an atom momentum 
p in the laboratory frame will assume a modified momentum 

Plascr = P - mvdex (20) 

in the laser frame. As a consequence, a wavevector kg in the laboratory frame will change to 

kojascr = PL/?i = ko - mvdex/h (21) 

in the laser frame. For the sake of comparison with our early results, we set kg = 0, and an initial internal state same 
as that in Eq. p^ . Thus the total wavefunction in the laser frame is given by 

l^iascr^, = J ^ e-'^^ (22) 

The state in Eq. (|22|) in the laser frame is not an eigenstate of the Hamiltonian. Nevertheless, it can be written as 
a superposition of two eigenstates |^'^'+) and I^E"^'^) defined previously by Eq. (fTTI) . Its time evolution can then be 
obtained by use of the two energy eigenvalues associated with \'^^'^) and \'i>^'^) (with k — kojascr)- Specifically, 

lascr/.NN _ 1 ( 1 I -i^^x -iE+t/h 



|^lascr(^)^, = - f ^ ) e"^ — 



2 I -1 



^-^^^^-^E-t/h^ ^23) 
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where and are the two eigenvalues [(/cq, laser i K)h]'^/2m, with fcojasor — "'^^x^' can obtain the same 

expression for the wavefunction in the laser frame more formally, by use of the well-known Galilean transformation 

for wavefunctions, namely, \^{t)) = |vI/'^>5<=''(t))e'~B~^+*TB~*. 

Interestingly, the two phase factors in Eq. (|23p. i.e., e~*^+*/'' and g"'^-*/'', are purely dynamical phase factors. 
This is in contrast to our early perspective in the laboratory frame (where a geometric phase may arise). The two 
group velocities corresponding to the two components of the state in Eq. (1^^ . both along the x direction, can then 
be obtained as follows, 

Vd + hK/m; 

Vd — hn/m. (24) 

We can now return to the laboratory frame. Because the laboratory frame is moving with the laser frame at a velocity 
—Vdfix, the two group velocities in the laboratory frame are also along the x direction, and they are given by 

^lab ^ v^^'"' + Vd^nK/m; 

^lab ^ v^^"""' + Vd = -hK/m. (25) 

This result, which is based entirely on dynamical phase consideration, is exactly the same as those in Eq. (jl9p 
obtained by considering the integral of the non-Abelian vector potential in the laboratory frame. Consistent with 
a previous finding that a geometric phase may not be Galilean invariant (isj . it is still intriguing to see from our 
concrete example that the effect of a non-Abelian vector potential (which leads to a gauge-independent geometric 
phase for cyclic processes) may be interpreted totally as that of a dynamical phase in a different frame of reference. 
That is, the quantum phases and their consequences can be traced back to either an integral of the vector potential 
in laboratory frame induced by laser beam displacement or to the dynamical phases in the laser frame (i.e., the time 
integral of the energy eigenvalues). In retrospect, this is not entirely surprising. The phase due to the vector potential 
is generated by active parameter manipulation associated with laser-field displacement. In a frame of reference in 
which laser fields are not moving, there is no such phase induced: observable effects (such as wavepacket splitting) 
must be understandable in terms of other evolution effect, which is purely the dynamical phase here. 



^laser 



dE+/h 

-laser 

laser dE_/h 



^lascr 



dkQ, 



laser 



B. Formation of a square or triangular lattice of cold-atom sub-wavepackets 

When the laser beams are displaced along the path A ~ B in Fig. 1, an initial wavepacket can split into two parts 
with group velocities ihn/m in the x direction. Here we first generalize this scheme by displacing the laser beams 
along the square shown in Fig. 1. 

The two wavepackets shown in Fig. 2(a) are associated with the two components shown in the last line of Eq. 
(|15p . with the only difference being that for wavepacket simulations a Gaussian profile is imposed on the initial state 
and hence each component does not extend to infinity. Because our plane-wave consideration is seen to describe the 
actual wavepacket dynamics very well, we continue to use plane waves to understand the process along the square. In 
addition, because the two wavepackets in Fig. 2(a) are already separated when the point B is reached, we only need 
to treat each term in Eq. ([15]) separately, rather than treating the two sub-wavepackets as a whole. 

Consider then the second stage of a clockwise laser beam displacement along the square, which is along the path 
B — C shown in Fig. 1, with a velocity Vd in the —z direction. The effect of the non-Abelian vector potential can be 
obtained from 

^-G.(:;), (26) 

where a 2 x 2 matrix Gz is defined. This effect now reflects the z-component of the non-Abelian vector potential. 
Because of the special form of the dark states that have two wavevectors k\. and k', here the parameter manipulation 
is always noncyclic, in the sense that the dark states do not return to their initial form after the laser displacement. 
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FIG. 2: (Color Online) Numerical simulation of the dynamics of tripod-scheme cold-atom wavepackets when all the three laser 
beams are slowly displaced along a square shown in Fig. 1. The calculations are based on the full Hamiltonian — ^^2m — ^ 
Hrwa,4 + Vb (In real units, we take m — lO"'^^ Kg, k ~ 10^ m~^). (a)— >-(d): wavepacket densities when the points B, C, D 
and A are reached, respectively, in the case of a clockwise displacement of laser beams. The time intervals for A — B and B — C 
are 20m/ (Tin^) ~ 0.03s and those for C — D and D — A are 4Qm/ {TiK^) ~ 0.06s. (e)— ^(h): wavepacket densities in the case of 
counter-clockwise laser beam displacement, when the points D, C, B and A are reached, respectively. The time intervals for 
displacement along A — D, D ~ C and C — B are 0.03s and that for D — A is 0.06s. The initial state is taken as that shown in 
Eq. (|12|) . X and z are in units of 1/k and a — 100/k. 



As such, the evolution matrix as a solution to Eq. (|26|) is necessarily representation-dependent. For this reason, in our 
fixed representation for the two dark states \Di) and 1 1)2), we do not call the evolution matrix determined by Eq. (pS)) 
as a non-Abelian geometric phase, but only as an effect due to the non-Abelian vector potential defined above. (To 
remove such a representation-dependence one can adopt a unique gauge proposed in Ref. ^19.] based on the concept of 
parallel frames via re-defining the dark states. However, it is found that the resulting non-Abelian gauge potential is 
too complicated. Alternatively, we can simply remove the factor e~"* ^ from the dark states defined previously. With 
the newly defined dark states, our parameter manipulation will be cyclic and the obtained solution to Eq. (|26p will 
have a desired gauge- invariance so long as the displacement along z is a multiple of 27r/fcJ.. But most importantly, 
irrespective of what representation we use, the final solution to Schrodinger's equation remains unchanged). 

We then proceed by assuming that the initial internal state is either the first or the second component of Eq. (jlSp . 
Apart from an overall phase when the point B is reached, the first component of Eq. will evolve to 



when point C is reached, where (pd is again a common dynamical phase for the internal states (1,0)^ and (0, 1)"^. 
To seek the group velocities for the two new components in Eq. (|27l) . we also need to find the derivatives of the 
two frequencies -iiVdH with respect to the wavevector ko in the laboratory frame. In the same manner as we derive 
Eqs. (jl8p and (jl9p . we may consider a ^fco in the z direction and take advantage of a frame moving at the velocity 
hSkoBz/m relative to the laboratory frame. Denote again Vd as the velocity of the laser beams in the frame with a 
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FIG. 3: (Color Online) (a) A schematic plot showing that the laser-beam displacement is along an equilateral triangle with the 
side length a = 100/fi:. The laser beams are moved either clockwise or counter-clockwise from the starting point A. The initial 
wavepacket is placed at the origin with the internal state taken the same as that shown in Eq. (|12p . (b) The final distribution 
of sub-wavepackets, obtained from numerical simulation of the tripod-scheme cold-atom dynamics when all the laser beams 
have been slowly displaced clockwise along the triangle for one cycle. The calculations are based on the full Hamiltonian 
~%m — HnwAA + Vs (In real units, we take m = 10~^^ Kg, k ~ 10^ m~^). The time interval to move across each side of the 
triangle is 20m/{hK,^) ~ 0.03s. x and z are in units of (c) Same as in panel (b), but the displacement of the laser beams 
is counter-clockwise. 



zero wavevector for the spatial part of the total wavefunction, we obtain 

(28) 



dvd h 
dko m ' 



and the two new group velocities along the z direction, 



lab ^ djvdK) dVd ^ 

dvd ' dko 
lab _ d{-VdK) did 



These two group velocities correspond to the internal states (1,0)^ and (0, 1)"^, respectively. It can then be predicted 
that the wavepacket associated with the first term in the last line of Eq. ([T5|) will further split into two parts in the 
z direction, due to the different group velocities ±hK/m. In the same manner we can predict that the wavepacket 
associated with the second term in the last line of Eq. (jTS]) will also split into two parts. As a consequence, when 
the point C is reached, there should be altogether four sub-wavepackets located ata; = ± — — , z = ± — — . The 
numerical results shown in Fig. 2(b) nicely confirm our predictions. 

As the laser beams are displaced further along the square shown in Fig. 1, in principle one can split the wavepacket 
into 2" copies on the x — z plane, forming a beautiful square lattice of matter waves. Figure 2(c) and 2(d) show 
the numerical results of the distribution of sub-wavepackets when the points D and A are reached, respectively, in 
a clockwise order. Remarkably, when the laser beams are moved back to the initial position A, there are altogether 
16 sub-wavepackets formed. It is interesting to examine the relative quantum phases between each sub- wavepacket. 
If their relative phases can be easily calculated and are independent of the details of the manipulation process, then 



9 



a coherent lattice of matter wavepackets is created and they should be useful for atom optics applications such as 
atomic interferometry. This is indeed the case here. In particular, because of a common dynamical phase (pd, the 
relative phases between the sub-wavepackets can be easily determined by examining the imprinted geometric phases, 
say zLvdKt shown in Eqs. ((T5|) and ((27)) . Take the four sub-wavepackets in Fig. 2(b) as an example. Other than a 
common overall phase, their quantum phases are found to be —2aK, 0, 0, and 2aK, for the upper-left, upper-right, 
lower-left, and lower-right sub-wavepackets, respectively. 

It should also be pointed out that, after displacing the laser beams along a square, a closed path in the parameter 
space is formed. As such, the overall effect due to the non-Abelian vector potential becomes that of a gauge- 
independent non-Abelian geometric phase. Further, because this manipulation is achieved by displacing the laser 
beams along different straight lines, our discussions above also indicate that this geometric phase effect can be equally 
understood as arising from the associated dynamical phases in several laser frames. 

The formation of a square lattice of cold-atom sub-wavepackets makes it clear that it should be equally possible 
to form other types of lattices. A simple consideration shows that if we move the laser beams along an arbitrary 
direction, the wavepacket splitting, if it occurs, must be along the same direction. As a result, if we move the laser 
beams along an equilateral triangle, it should be possible to form a triangular lattice of cold-atom sub-wavepackets. 
This is confirmed by our results shown in Fig. 3. In particular. Fig. 3(b) [Fig. 3(c)] depicts the distribution of 
sub-wavepackets, with unequal weights, after the laser beams are moved one clockwise (counter-clockwise) cycle along 
an equilateral triangle shown in Fig. 3(a). By repeating this process more sub-wavepackets can be generated. 



C. Non-Abelian Aharonov-Bohm effect 

In Refs. [13, [2l|, the non-Abelian Aharonov-Bohm effect is demonstrated by showing different internal states 
corresponding to two different paths from a common starting point to a common ending point, or by showing the 
non-commutability between different navigation paths. To make connection between our results and the context of 
non-Abelian Aharonov-Bohm effect, here we consider in detail what happens to a wavepacket if the laser beams are 
displaced counter-clockwise along a square, with the same initial state in Eq. (jl2l) . In particular, along the path 
A — D (see Fig. 1) the laser beams are displaced in the —z direction. The internal state will then evolve in accord to 
Eq. (|26)) plus a dynamical phase (j)d- When the point D is reached, then the system is in the state 

l^Dito)) = e-''^^^^'" (1)^'^" (30) 



Following the same procedure for deriving the group velocities in Eq. (j29[) . we obtain that the phase factor —VdKtD = 
—Ka will induce one single group velocity hn/m in the z direction. So when the point D is reached the wavepacket 
does not split during the passage A — D but moves with a velocity fin/m in the z direction. The numerical result 
shown in Fig. 2(e) confirms this analytical result. 

The next stage of a counter-clockwise laser beam displacement is along the path D — C shown in Fig. 1. Apart 
from an overall phase, the total wavefunction when the point C is reached can be obtained from the phase induced 
by the non-Abelian vector potential and a dynamical phase previously defined, i.e., 

\^c{tc)) = e'^^vd{tc-tn) (^^^^ ^^^d (31) 
1 / 1 A ^iVdK,{tc-to)^itl>d ^ i I I g-iVdK.(tc-tD) gi4>d 



2 V 1 / 2 V -1 

Associated with the two internal states "^(l, 1)"^ and -^(1, ~l)"^i we find that the two group velocities are zLTiK/m, 
both in the x direction. Thus an initial wavepacket will split into two parts in the x direction when the point C 
is reached in a counter-clockwise order. Numerical result in Fig. 2(f) again confirms this. Further splitting of the 
wavepackets is shown in Fig. 2(g) and 2(h). 

By comparing the results associated with clockwise laser beam displacement and those associated with counter- 
clockwise displacement, we can now clearly see a non-Abelian Aharonov-Bohm effect. In particular, the wavepacket 
splitting behavior due to displacement along the path A — B — C (four equal- weight copies in Fig. 2(b)) is different 
from that along the path A — D — C (two equal- weight copies in Fig. 2(f)). This dramatic difference shows that a 
non-Abelian Aharonov-Bohm effect can manifest clearly on the translational motion of cold atoms, rather than on 
their internal states considered in Refs. [20l. [2li. 
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FIG. 4: The scheme of laser beam displacement along a circle with a radius tl, with a starting point E. The atom wavepacket 
will be forced to dance along a circle with a radius r,4 during the process of circular laser beam displacement. The initial group 
velocity of a cold-atom wavepacket in the z direction will be changed to v/ when the point F is reached. 



IV. RELATION BETWEEN A NON-ABELIAN GEOMETRIC PHASE AND AN ABELIAN 

GEOMETRIC PHASE 

Results in the previous section show that by considering different frames of reference, new insights into the dynamics 
induced by laser beam displacement may be obtained. In this section we extend this idea and attempt to make a 
connection between the non-Abelian vector potential induced phase and an Abelian geometric phase. As a concrete 
example, we consider a situation where all the laser beams are slowly moved along a circle of a radius ■ After one 
cycle, a closed path in the parameter space is formed, and as such the integral of the non-Abelian vector potential 
gives rise to a non-Abelian geometric phase that is gauge invariant. 



A. Perspective from the laboratory frame 

As shown in Fig. 4, we slowly displace the laser beams from the starting point E along the circle Cl clockwise in 
the laboratory frame, with the initial {t — 0) direction of the laser displacement in the z direction. The moving speed 
of the laser beams is still denoted by Vd- The initial state in the laboratory frame is again chosen as that in Eq. (|12p . 

In a circular laser displacement process, the direction of the laser displacement always changes, with an angular 
velocity Vd/rL. At time t, the laser beams are moving in the direction sin(^i)ea; -I- cos(^t)e2. Let s — Vdt, at time t 
we have 

dx . Vd,. 

— = sm{—t), 
as rL 

— = cosi — t). (32) 
as rL 

Combining Eqs. ([13]), (|26p and (15^ . one can easily find the evolution of (ci, C2) at time t satisfies 

As explained in the previous section for kp — 0, the dynamical phase 0^ given by Eq. (|14l) is just an overall phase for 
any internal state in the dark subspace. 4>d can hence be neglected. 

Interestingly, though Eq. (p3)) clearly describes the effect of a changing non-Abelian vector potential, it however 
assumes the very same form as a time-dependent Schrodinger equation of a two-level system. As a result its solution 
can be obtained by treating Eq. ([55]) as a Schrodinger equation for a time-dependent pseudo Hamiltonian 

/ cos(^i) sin(^t) \ , , 

i^.W = -ta ( ^^g,; -cos(|^t)j- (34) 
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The slow movement of the laser beams induces a slow change of this pseudo Hamiltonian. The eigenstates of Hp{t) , 
denoted |1(2)), are 

|2) - |.gk(t)>, (35) 
where \g^l^'^ ^) are the internal state defined in Eq. (jlip. with the angle between the x axis and k(t) given by 



The energy eigenvalues of |1(2)) as eigenstates of Hpit) are 

El = -hnvd; 

E2 = TiKVd- (37) 

Since Ei — E2 ^ 0, the two eigenstates |1(2)) are necessarily non-degenerate. Further, because at i = 0, the initial 
state we adopt here (given by Eq. (fT2)) ) is exactly the eigenstate |1) with 1^915.(4) = 7r/2, the time evolution of this 
initial state with a slowly changing pseudo Hamiltonian Hp{t) becomes a non-degenerate adiabatic problem. That is, 
so long as the change of this pseudo Hamiltonian is sufficiently slow, then the adiabatic theorem for a non-degenerate 
spectrum guarantees that the time-evolving state will remain on the instantaneous eigenstate l^^^j^) afterwards. 
Roughly speaking, a sufficiently slowly changing "iJp(i)" requires [l^ 

«1. (38) 



\Eiit) - E2{t)\ 

Substituting Eqs. (|35|) . (jlip and (|37|) into Eq. (|38|) . we find that the above adiabatic condition, when applied to the 
pseudo Hamiltonian H.p{t), reduces to 

rL»-. (39) 

K 

If this additional condition is satisfied in the laser field manipulation, then according to the standard adiabatic theorem 
for a non-degenerate spectrum, at time tp when an arbitrary point F along the circle is reached, the initial state 
IS'k(t=o) ('''k(t=o) — §)) — (IjO)"^ will evolve to the instantaneous eigenstate |.gi^(tF)) multiplied by a geometric phase 
factor, namely, 

where a is the angle between x axis and the tangent line of the circle at the ending point F shown in Fig. 4. The 
phase j3 is given by the following expression that is characteristic of an Abelian geometric phase, 

P = J^~j9i\^\9i)-d^. (41) 

Because this geometric phase arises from considering a pseudo Hamiltonian, it could be called as a pseudo Berry-like 
Abelian geometric phase. Recalling that the two-component final state given by Eq. PO]) after moving the laser 
beams along one cycle is the result of non-Abelian geometric phases ([33|) imprinted onto the initial state (1,0)^, all 
the components 1 — ie^'^^, —i + e*''''' and the phase /3 comprise the total effect of the non-Abelian vector potential. In 
this sense /3 is just one part of a total non-Abelian geometric phase. 

Figure 4 illustrates what can be predicted from this picture afforded by the pseudo Hamiltonian defined above. 
When t — the laser beams are displaced in the z direction and the initial state is 15^(4— o)('^k(t=o) — f ))i when the 
ending point F is reached at t = tp, the laser beams are displaced in the indicated a direction and the final state 

\9k{t=tF)('fiKt=tF] ~ '^))- '^he group velocity of the time-evolving state can be obtained using the same method 
as that used for deriving Eqs. (IT5|) and (I^Hl) . It is found that the magnitude of the group velocity is fixed at Tin/ra, 
but its direction changes from 'y9k(t=o) — '''/2 to ^P\^(t=tF) ~ Thus, the instantaneous laser displacement direction 



12 




-50 50 -50 50 
X X 

FIG. 5: (Color Online) Numerical simulation of the wavepacket dynamics of tripod-scheme cold atoms, with the involved laser 
beams being slowly displaced along a circle shown in Fig. 4. The white circle illustrates a possible circular trajectory traced 
by the center of the moving wavepacket. x and z are in units of See the text for details. 

always coincides with the group velocity of the wavepacket. Now if the speed of laser displacement Vd is a constant, 
then one deduces that the atom wavepacket must dance on a circle Ca with a different radius r^i, which satisfies 

rA ^ hn/m 

This also suggests that one can easily change the final group velocity of a wavepacket. For example, if we need to 
change the initial group velocity in the z direction to the a direction as shown in Fig. 4, we can displace the laser 
beams along the circle Cl to point F and then along the tangent line at point F. 

All these predictions have been verified by our numerical simulations. In particular, in our simulations the initial 
state is chosen as a Gaussian wavepacket instead of a plane wave considered in Eq. (|12p . As an example we initially 
locate the wavepacket a.tx — — 50/k, z = 0. The laser beams are displaced with a velocity l.bfiK/m along the circle Cl 
of radius vl « 75/ k. As shown in Figs. 5(a) and 5(c), the wavepacket indeed dances on a circle of radius va » 50/ k, 
satisfying the relation By contrast. Fig. 5(b) depicts the wavepacket a,t t = 50to/(?ik^), if the laser beams are 

displaced not along the circle Cl but in the z direction. The vertical displacement of the wavepacket indeed shows 
that the initial group velocity of the initial state is in the z direction. Figure 5(d) depicts the wavepacket after the 
lasers are displaced along the circle Cl to point F and then along the tangent line for a duration of 50m/(?iK^). The 
result in Fig. 5(d) demonstrates that the final group velocity is indeed in the direction of the tangent line at point F. 



B. Perspective from a different frame of reference 



In the previous subsection we have argued that one part of a non-Abelian geometric phase, namely, the /3 phase 
factor in Eq. I|4ip . can be regarded as a pseudo Abelian geometric phase. In this section, we revisit this connection 
between a non-Abelian geometric and an Abelian geometric phase in the laser frame of reference. As shown below, 
in this new frame of reference the pseudo Abelian geometric phase (3 becomes a truly Abelian geometric phase, i.e., 
an Abelian geometric phase without introducing a pseudo Hamiltonian. 

As seen before, in the laser frame the wavevector for the spatial part of the total wavefunction becomes kojasor = 
—Vdm/h. Note that now the direction of the laser beam displacement is continuously changing with time. That is. 



Vd = Vd 



cos( — t)ez + sm( — t)ex 

TL TL 



(43) 
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with its magnitude fixed at Vd- Note also that in the laser frame the wavevector for the translational part of the total 
wavefunction should be given by ko.iasci - One then finds that in the laser frame the total Hamiltonian in Eq. (jlOp 
can be explicitly written as 

.fF ^^(^0 laser + '«^) fcosi^t) sin(^i) \ , , 

The above Hamiltonian assumes the very same form, except the constant term ?i^(fco laser + K^)/2m, as the pseudo 
Hamiltonian in Eq. (|34|) . Clearly then, the adiabatic evolution of the eigenstate of H]^ will require exactly the same 
adiabatic condition as Eq. (|39|) . Under this adiabatic condition, the final result of the adiabatic process viewed in 
the laser frame should be given by the state in Eq. (j40p multiplied by a component due to ko.iascr, i-C, 

l*/,iascr> = I ( ) e'^'e^'^o.'--^. (45) 

Returning to the laboratory frame, the wavevector for the spatial part of the total wavefunction changes from ko, laser to 
kg = 0, so the above state |5'/,iascr) reduces exactly to that obtained in Eq. (|40)) . The two perspectives from different 
frames of reference are thus in agreement. Because here stands for a real Hamiltonian for the translational 
motion in the laser frame, the pseudo geometric phase /3, which is part of the total non-Abelian geometric phase, 
can now be understood as a truly Abelian geometric phase in the laser frame. Thus an interesting relation between 
non-Abelian and Abelian geometric phases is established, i.e., in some cases, part of non-Abelian geometric phase can 
be interpreted as an Abelian geometric phase from a different frame of reference. 



V. CONCLUDING REMARKS 

To conclude, we have considered two interesting extensions of our previous study [Tsj of quantum control of matter 
waves in tripod-scheme cold-atom systems. The results have enhanced our understandings on a number of issues and 
should motivate further experimental and theoretical interests in the dynamics of tripod-scheme cold atoms. 

In particular, by displacing the laser beams along a square, it is shown that the effect of a non-Abelian vector 
potential can split an initial wavepacket into 2" copies that form a coherent square lattice. It is also shown that 
even a triangular lattice of cold-atom sub-wavepackets can be formed by moving the laser beams along a triangle. 
Coherence between the sub-wavepackets thus generated makes them potentially useful for atom optics applications. 
We further discussed how to clearly manifest a non-Abelian Aharonov-Bohm effect by comparing the wavepacket 
profiles obtained by clockwise and counter-clockwise paths along a square. Examining the same process from a 
different frame of reference that moves with the laser beams, we find a remarkable connection between an integral of a 
non-Abelian vector potential (which becomes a gauge-independent non-Abelian geometric phase for cyclic processes) 
and simple dynamical phase expressions in the laser frame. 

By displacing the laser beams along a circle, we show that a non-Abelian vector potential effect can force the 
atom wavepacket to dance on a circle of a different radius and can also redirect the propagation direction of that 
wavepacket. Comparing the perspective in the laboratory frame and in the frame that moves with the laser beams, 
it is shown that one part of the associated non-Abelian geometric phase in one frame of reference can be understood 
as an Abelian geometric phase in another frame of reference. 
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